Overlap Dirac operator at nonzero chemical potential and random matrix theory 
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We show how to introduce a quark chemical potential in the overlap Dirac operator. The resulting 
operator satisfies a Ginsparg- Wilson relation and has exact zero modes. It is no longer 75-hermitian, 
but its nonreal eigenvalues still occur in pairs. We compute the spectral density of the operator 
on the lattice and show that, for small eigenvalues, the data agree with analytical predictions of 
nonhermitian chiral random matrix theory for both trivial and nontrivial topology. 
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Recent years have seen great advances in two areas that 
at first sight seem to be totally unrelated: (i) the study of 
nonhermitian operators in the natural sciences l,] and (ii) 
the problem of chiral symmetry on the lattice j2| ■ In this 
article we focus on a problem in which both of these areas 
are relevant, namely quantum chromodynamics (QCD) 
at nonzero baryon (or quark) density, which is important 
for the study of relativistic heavy- ion collisions, neutron 
stars, and the early universe [2|. 

If a quark chemical potential ^ is added to the QCD 
Dirac operator, the operator loses its hermiticity prop- 
erties and its spectrum moves into the complex plane. 
This causes a variety of problems, both analytically and 
numerically. Lattice simulations are the main source of 
nonperturbative information about QCD, but at fi ^ 
they cannot be performed by standard importance sam- 
pling methods because the measure of the Feynman path 
integral, which includes the complex fermion determi- 
nant, is no longer positive definite. While a generic solu- 
tion to this so-called sign problem is unlikely to be found 
|j|, a number of recent works have been able to make 
progress by circumventing the problem in various ways 
la, lo, I2I • These methods all agree on the transition tem- 
perature from the hadronic to the quark-gluon phase in 
the regime /x/T < 1 0. 

A better analytical understanding of QCD at very high 
baryon density has been obtained by a number of meth- 
ods [S|, and the QCD phase diagram has been studied 
in model calculations based on symmetries J9|. Chiral 
random matrix theory (RMT) [lOJ, which makes exact 
analytical predictions for the correlations of the small 
Dirac eigenvalues, has been extended to /i 7^ [ll|, and 
a mechanism was identified [l3| by which the chiral con- 
densate at ^ 7^ is built up from the spectral density of 
the Dirac operator in an extended region of the complex 
plane, in stark contrast to the Banks-Casher mechanism 
at /i = 0. 

A first comparison of lattice data with RMT predic- 
tions at /i 7^ was made in Ref. [13 using staggered 
fermions. One issue with staggered fermions is that the 
topology of the gauge field is only visible in the Dirac 
spectrum if the lattice spacing is small and various im- 
provement and/or smearing schemes are applied [lj| . To 



avoid these issues, we would like to work with a Dirac 
operator that implements a lattice version of chiral sym- 
metry and has exact zero modes at finite lattice spacing. 
The overlap operator [l3| satisfies these requirements at 
/i = 0. In the following, we show how the overlap opera- 
tor can be modified to include a nonzero quark chemical 
potential [la, Uj • We then study the spectral proper- 
ties of this operator as a function of fi and compare data 
from lattice simulations with RMT predictions. As we 
shall see, the overlap operator has exact zero modes also 
at nonzero /x, which allows us, for the first time, to test 
predictions of nonhermitian RMT for nontrivial topology. 
We begin with the well-known definition of the Wilson 
Dirac operator Dw including a chemical potential /k [l^ , 



Dwip) 



nY,{Tt + T-)~K{e^T+ + e-^T-, 



i=l 



{T'^)yX = (1 ± lu)U±y{x)5y^x±0 , 



(1) 



where k — l/{2m\y + 8) with the Wilson mass mw, the 
U e SU(3) are the lattice gauge fields, and the 7^ are 
the usual euclidean Dirac matrices. Unless displayed ex- 
plicitly, the lattice spacing a is set to unity. The overlap 
operator is defined at /x = by [l5l | 



£'ov(0) 



75 £(75-0^/(0)) 



(2) 



where e is the matrix sign function and 75 = 71727374- 
niiY must be in the range (—2,0) for DqvIO) to describe 
a single Dirac fermion in the continuum. The properties 
of -Dov(O) have been studied in great detail in the past 
years. In particular, its eigenvalues are on a circle in 
the complex plane with center at (1,0) and radius 1, its 
nonreal eigenvalues come in complex conjugate pairs, and 
it can have exact zero modes without fine-tuning. Z3ov(0) 
satisfies a Ginsparg- Wilson relation [l9j of the form 



{D,-f5}^Djr,D^ 



(3) 



We now extend the definition of the overlap operator to 
/i 7^ 0. The operator Dw{0) m Eq. (0) is 75-hermitian, 
i.e., 75i3vy(0)75 = D{y(0), and therefore the operator 
"fbDwiO) in the matrix sign function is hermitian. How- 



ever, for n j^ 0, DwifJ-) is no longer 75-hermitian. Defin- 
ing the overlap operator at nonzero /i by 



DovifJ-) = l + "f5£{l5DwitJ-)) 



(4) 



we now need the sign function of a nonhermitian matrix. 
In general, a function / of a nonhermitian matrix A can 
be defined by a contour integral. A more convenient ex- 
pression can be obtained if A is diagonalizable. In this 
case we can write A = UAU^^, where U G G\{N, C) with 
N = dim(A) and A = diag(Ai, . . . , Aat) with A^ € C. 
Then f{A) = C//(A)C/"\ where /(A) is a diagonal ma- 
trix with elements /(Ai). The sign function can be de- 
fined by [13 



s{A) = Us[gn{ReA)U- 



(5) 



This definition ensures that e'^{A) = 1 and gives the 
correct result if A is real. An equivalent definition is 
e{A) = A{£-Y^l'^ EJ]. Eqs. gj) and 10 constitute our 
definition of _Dov(m)- The sign function is ill-defined if 
one of the \i lies on the imaginary axis. Also, it could 
happen that ^^Dm(^[i) is not diagonalizable (one would 
then resort to a Jordan block decomposition). Both of 
these cases are only realized if one or more parameters 
are fine-tuned. This is unlikely to happen in realistic 
lattice simulations, and we therefore ignore these issues. 
It is relatively straightforward to derive the following 
properties of I?ov(/^)- 

• Dov{fi) is no longer 75-hermitian because for a non- 
hermitian matrix A we generically have e^(A) 7^ £{A). 
Instead, DovifJ^) satisfies j5Dov{fJ-)j5 = DI^{~^). 

• Z?ov(/^) still satisfies the Ginsparg- Wilson relation (O 
because of e'^{A) — 1. Thus, we still have a lattice 
version of chiral symmetry, and the operator has exact 
zero modes without fine-tuning. 

• The eigenvalues of Z3ov(m) that are not equal to or 2 
no longer come in complex conjugate pairs, but every 
such eigenvalue A (with eigenfunction ip) is accompa- 
nied by a partner A/(A — 1) (with eigenfunction jstp)- 

• At /Lt = 0, the mapping X ^ z = 2A/(2 — A) projects 
the eigenvalues from the circle to the imaginary axis. 
At /i 7^ 0, the same mapping projects the eigenvalues A 
and A/(A — 1) to a pair ±z, which is the same pairing 
as for /i / in the continuum. 

• The eigenfunctions of Dov(a') corresponding to eigen- 
value or 2 can be arranged to have definite chirality. 
The numbers nf of modes with A = 0, 2 and (75) = ±1 
satisfy Uq — Uq ~ — (nj — f^^ ) . (Without fine-tuning, 
we even have n^ = n^ = or tiq — n^ = 0, which is 
always the case in our lattice data below.) 

• Generically, _Dov(m) is not normal, i.e., DD'^ 7^ D'^ D. 

We now turn to our (quenched) lattice simulations. 
The computation of the sign function of a nonhermitian 
matrix is very demanding. We are currently investigat- 
ing various approximation schemes [22l |. but in this ini- 
tial study we decided to compute the sign function and 
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FIG. 1: Spectrum of Dov(m) for /i = and fj, = 0.3 for a typ- 
ical configuration. The figure on the right is a magnification 
of the region near zero. 



to diagonafize Dov(a') exactly using LAPACK. For the 
comparison with RMT we need high statistics, which re- 
stricts us to a very small lattice size. We have chosen 
the same parameter set as in Ref. [23 to be able to com- 
pare with previous results at ^ = 0. The lattice size is 
V — 4"^, the coupling in the standard Wilson action is 
(3 — 5.1, the Wilson mass is mw = ^ 2, and the quark 
mass is mq — 0. The small lattice size forces us to use 
such a strong coupling to stay in the ergodic regime of 
QCD, see Eq. (fTT|) below, but we expect the conclusions 
of this paper to remain valid at weaker couplings and cor- 
respondingly larger lattice sizes. The choice of ^m^r — —2 
is motivated by concerns about the locality of Dov(m)i al- 
though at y = 4^ this question is largely academic. 

Our parameters are given in Tabled In Fig.^wc show 
the spectrum of Dov(m) foi' a typical configuration for fj, = 
and n = 0.3. As expected, we see that the eigenvalues 
move away from the circle as fi is turned on. Another 
observation is that the number of zero modes of Dov (m) 
for a given configuration can change as a function of /z, 
see Table 13 This can be understood from the relation 
between the anomaly and the index of Z?ov |2J, |25j , 



'tr(75£)ov) = 2index{Dov) , 



(6) 



which we can show to remain valid at /i 7^ 0. Using 
tr(75£'ov) = tr[75 + e{j5Dw)] = ti[e{'y5Dw)] and the 
fact that the eigenvalues of the sign function are -1-1 or 



- 1 , one has n 



w 



2 index(£'o 



where n^ denotes 



the number of eigenvalues oij^Dw^fx) with real part ^ 0. 
As ^ changes, an eigenvalue ofj^Dw can move across the 



imaginary axis. As a result. 



,w 



changes by 2, and 



TABLE I: Number of configurations and distribution of the 
number f of zero modes of Dov(m) for various values of /i. 



fia 



# config. 



P{i> = 0) P{v = 1) Pif = 2) 
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0.476 


0.465 
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FIG. 2: Density of the small eigenvalues of Dov{y) in the complex plane (after projecting A — > 2A/(2 — A) and rescaling by VTj) 
for (from left to right) ^ — 0.1,0.2,0.3, 1.0. The histograms are lattice data for v — and u = 1, and the solid lines are the 
corresponding RMT prediction of Eq. 0, integrated over the bin size. Top: cut along the imaginary a^cis, middle: cut along 
the real axis, bottom: cut parallel to the real axis. The vertical lines indicate the fit interval. For fi = 1.0, the eigenvalues are 
rescaled by VT,/2^/a, the fit is a one-parameter fit to Eq. po^. and in the bottom plot the data are integrated over the phase. 



thus index(Dov) changes by 1, which explains the obser- 
vation. We believe that this is a lattice artefact which 
will disappear in the continuum limit. (Also, the index 
theorem is violated at /3 — 5.1, but for the comparison 
with RMT only the number of zero modes matters.) 

The spectral density of DqvIm) is given by p°'^{Xr, Xi) = 
{J2k ^(-^ ~ ^k)) with X = Xr + iXi, where the average is 
over configurations. The claim is that the distribution of 
the small eigenvalues of Dov{fi) is universal and given by 
RMT. The RMT model for the Dirac operator is [ul 



D 



RMT 



(a^) = 







H 



^^ 



(7) 



where VK is a complex matrix of dimension n x (n + v) 
with no further symmetries (we take v > Q without loss 
of generality). The matrix in Eq. (|7|) has v eigenvalues 
equal to zero. The spectral correlations of Drmt(m) on 
the scale of the mean level spacing were computed in 
Refs. [2^ I23, 1231 . In the quenched approximation, the 
result for the microscopic spectral density reads 



„RMT 



{x,y) 



2Tra 

r-l 



e ^°' 



-K, 



x^ + j/^ 



4a 



tdtt 



-lat^ 



luitzt 



(8) 







where z = x + iy, I and K are modified Bessel func- 
tions, and a = ^^ f^V . To compare the lattice data to 



this result, the eigenvalues Xu need to be rescaled by a 
parameter \/VYi which is proportional to the mean level 
spacing near zero. 
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(9) 



S and /tt are low-energy constants that can be obtained 
from a two-parameter fit of the lattice data, Eq. Q, to 
the RMT prediction, Eq. JSJ. (The normalization is fixed 
by J dx dy p°^ (x , y) = 12^ and does not introduce an- 
other parameter.) For x <$^ a, Eq. (|Hl simplifies to |22J 



P. 



RMT 



(x^y) 



27ra 



K, (C) /. (0 with e 



\z\2_ 

4a 



(10) 



A fit of Eq. Q to Eq. ifTUIl then only involves the single 
parameter T^/ f^- 

In Fig. [21 we compare our lattice data to the RMT pre- 
diction. We display various cuts of the eigenvalue density 
in the complex plane as explained in the figure captions. 
The data agree with the RMT predictions within our 
statistics. S and /^ were obtained by a combined fit to 
the z^ = and v = \ data for all three cuts and are 
displayed in Tabic UTI (These numbers have no physical 
significance at (3 = 5.1.) Separate fits to the 1/ = and 
V = 1 data give results consistent with those in Table ITll 
For /x = 1.0 we have fitted to Eq. H10|) since the distribu- 
tion of the small eigenvalues is radially symmetric up to 



TABLE II: Fit results for E and /^ (see text). 
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0.67 




0.2 
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0.320(4) 
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1.03 


|3 


1.0 


- 


- 


0.603(18) 


0.42 


[4 



1^1 ~ 0.7, and therefore S and /tt cannot be determined 
independently by a fit to Eq. ©. We note in passing 
that unfolding along the imaginary axis results in better 
agreement of the lattice data with RMT for larger values 
of y, but we shall not discuss this issue here since it is 
a finite-volume effect that will be unimportant on larger 
lattices. 

The RMT results are valid in the "ergodic regime" 



1 

1j 



(11) 



where L'^ = V, ttIt^ is the Goldstone boson mass, and 
A is the mass scale of the lightest non- Goldstone parti- 
cle. The rescaled eigenvalues z = XVS should therefore 
be described by RMT for |z| ^ /^V^ (the dimension- 
less Thouless energy). For our data, f^vV « 1 (this 
explains the choice of /3 = 5.1; for larger /3 this num- 
ber would be even smaller). Nevertheless, the data seem 
to be described by RMT quite well even a bit beyond 
this expectation. Also, /i = 1.0 corresponds to fiL — 4, 
which violates the inequality Hll|l . The fact that RMT 
still works below the Thouless energy in this case means 
that the zero-momentum modes decouple from the par- 
tition function. A detailed study of the range of validity 
of RMT at /i 7^ as a function of the lattice parameters 
will be the subject of further work (see also Ref. JO])- 

In summary, we have shown how to include a quark 
chemical potential in the overlap operator. The opera- 
tor still satisfies a Ginsparg- Wilson relation and has ex- 
act zero modes. The distribution of its small eigenvalues 
agrees with predictions of nonhermitian RMT for trivial 
and nontrivial topology. Our initial lattice study should 
be extended to weaker coupling, larger lattices, and bet- 
ter statistics. Work on approximative methods to enable 
such studies is in progress. For small volumes, reweight- 
ing with the fermion determinant should allow us to test 
RMT predictions for the unquenched theory 31j. 
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